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Abstract 

Recently, remarkable advances have been made in coupling a number of high-Q modes of nano- 
mechanical systems to high-hnesse optical cavities, with the goal of reaching regimes where quantum 
behavior can be observed and leveraged toward new applications. To reach this regime, the coupling 
between these systems and their thermal environments must be minimized. Here we propose a novel 
approach to this problem, in which optically levitating a nano-mechanical system can greatly reduce 
its thermal contact, while simultaneously eliminating dissipation arising from clamping. Through 
the long coherence times allowed, this approach potentially opens the door to ground-state cooling 
and coherent manipulation of a single mesoscopic mechanical system or entanglement generation 
between spatially separate systems, even in room temperature environments. As an example, we 
show that these goals should be achievable when the mechanical mode consists of the center-of-mass 
motion of a levitated nanosphere. 
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One of the most intriguing questions associated with quantum theory is whether effects 
such as quantum coherence and entanglement can be observed at mesoscopic or macro- 
scopic scales. As a first step towards resolving this question, recently much effort has been 
directed toward quantum state preparation of a number of high-Q modes of nano- and 
micro-mechanical systems - in particular, cooling such systems to their quantum ground 
state Reaching a regime where the quantum behavior of mechanical devices emerges 
is not only of fundamental interest, but could lead to new applications in fields such as 

n i 

ultra-sensitive detection [2| and quantum information science [3J. To reach this regime, it 
is critical that the thermalization and decoherence rates of these systems be minimized, 
by reducing the coupling to their thermal reservoirs. Thus far, this has necessitated the 
use of cryogenic operating environments. From an engineering standpoint, it would also be 
desirable to reduce the dissipation and thermalization rates of these systems through their 
clamping and materia, supports Q. so that these rates might approach then fundamenta. 
material limits [5fl. 

Here we propose a novel approach toward this problem, where the material supports are 
completely eliminated by optically levitating [f| a nano-mechanical system inside a Fabry- 
Perot optical cavity. Indeed, since the pioneering work of Ashkin on optical trapping of 
dielectric particles |f| (in the classical domain), it has been realized that levitation under 
good vacuum conditions can lead to extremely low mechanical damping rates 0, Q]. We 
show, however, that such an approach should also facilitate the emergence of quantum 
behavior even in room-temperature environments, when the particles are of sub-wavelength 
scale such that the effects of optical scattering become negligible. As a specific example, 
we show that the center-of-mass (CM) motion of a levitated nanosphere can be optically 



self-cooled 
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ll| to the ground state starting from room temperature. This system 



constitutes an extreme example of environmental isolation, as the CM motion is naturally 
decoupled from the internal degrees of freedom in addition to the external isolation provided 
by levitation. The long coherence time also allows for the preparation of more exotic states 
through coherent quantum evolution. Here, we consider in detail two examples. First, 



we describe a technique to prepare a squeezed motional state, which can subsequent 
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y be 
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mapped onto light leaving the cavity using quantum state transfer protocols 
Under realistic conditions, the output light exhibits up to ~15 dB of squeezing relative 
to vacuum noise levels, potentially making this system a viable alternative to traditional 
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techniques using nonlinear crystals 
shared between two modes of light 



16 1 . Second, we show that entanglement originally 
171 ] can be efficiently transferred onto the motion of 



two spheres trapped in spatially separate cavities, creating well-known Einstein-Podolsky- 
Rosen (EPR) correlations [18| between the mechanical systems. Ou r ap proach of optical 
levitation mirrors many successful efforts to cool jl^i I20I]. manipulate 21] and entangle 22] 
the motion of atoms and ions in room-temperature environments. At the same time, our 
system has a number of potential advantages, in that it enables direct imaging via strong 
fluorescence, exhibits large trap depths, and has a relatively large mass. We also note 
recent related experiments involving opto-mechanical "fluids" (with a continuous excitation 



spectrum rather than discrete modes) in the form of trapped, ultracold atomic gases [23|, |24|. 

Beyond the examples presented here, the use of a levitated device as an opto-mechanical 
system could open the door to many interesting opportunities. For instance, it should allow 
mechanical damping to approach fundamental material limits, potentially enabling the ex- 
ploration of nanoscale material properties. By levitating systems with internal vibrational 
modes, multiple modes could be optically addressed and cooled. In addition, the CM oscil- 
lation frequency can be tuned through the trapping intensity, allowing for adiabatic state 
transfer |25| with other modes or matching spatially separate systems for optical linking and 



entanglement generation [26|. Furthermore, this paradigm integrates nano-mechanics with 
many techniques for atomic trapping and manipulation, which can be further extended by 
levitating systems containing an internal electronic transition (e.g., a color center within a 
nano-crystal [stJ). Finally, as illustrated by squeezed light generation, engineering mechani- 
cal nonlinearities in conjunction with quantum state transfer yields a novel means to realize 
nonlinear optical processes. 

I. OPTICAL FORCES AND NOISE ACTING ON A DIELECTRIC SPHERE 

To illustrate our idea, we consider a sub-wavelength dielectric sphere interacting with 
two standing- wave optical modes of a Fabry- Perot cavity (Fig. [Ih.). One resonantly driven 
mode provides an optical dipole trap for the sphere. The second mode is driven by a weaker 
"cooling" beam, assumed to have a non-zero intensity gradient at the trap center, which 
provides a radiation pressure cooling force , lll[ . We discuss the cooling mechanism in 



the next section, while here we focus on the trapping potential and the noise forces acting 
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on the sphere. 

The trapping beam provides a gradient force similar to that used to "optically tweeze" 
small dielectric particles Q]. Considering a sphere whose radius is much smaller than the 
optical wavelength, r<^cA, its optical response is like that of a point dipole with induced 
dipole moment p in d = a in dE(x) and optical potential U opt (x) = — (l/4)(Re a in d)\E(x)\ 2 (see 
Appendix). Here x is the CM position of the sphere, a; n d = 3eoV^ (§r|) is its polarizabil- 
ity, V is the sphere volume, and e is the electric permittivity. Taking a standing wave 
E(x) = E cos kx (k=2ix/X), to lowest order near an anti-node the potential corresponds to 
a harmonic oscillator with mechanical frequency 

{Qk 2 I e — 1 \ 1//2 

uJ m = ° Re ^^ > (!) 

\ pc e + 2J 

where Jo is the field intensity and p is the mass density of the sphere. The total trap depth 
is Uq = (3/ ^/c)Re|=Tj. Typical trap depths and oscillation frequencies for a high-index 
material are plotted in Figs. [Tb,d. Frequencies of several MHz are achievable using 

an intra-cavity intensity of Io~l W//im 2 . The imaginary component of e characterizes 
optical absorption, which contributes to internal heating. For a material with ~10 dB/km 
propagation losses in bulk, intensities of Jo~10 W / pm 2 can be sustained without melting 
the sphere, due to blackbody re-radiation of the absorbed energy (see Appendix). With this 
in mind, we assume e is real in following discussions. 

The dominant noise forces acting on the sphere are collisions with a background gas 
and momentum recoil kicks due to scattered photons. In the Appendix, we show that the 
contributions from shot noise, blackbody radiation, and sphere anisotropy are negligible. 
Furthermore, the CM is de-coupled from the internal degrees of freedom and the sphere ef- 
fectively has no internal structure (as opposed to molecules, where the internal configuration 



can affect cooling efficiency [28| ) . In the regime where the molecular mean free path exceeds 
r, the background gas leads to a mean damping force dp/dt = — 7 9 p/2 with damping rate 
7 5 /2 = (8/ir)(P/vrp), where P, v are the background gas pressure and mean speed, respec- 



tively 221]. The random nature of the collisions also thermalizes the motional energy, at a 
rate given through the fluctuation-dissipation theorem by dE/dt = —j g (E — ksT), where T 
is the gas temperature. In particular, the characteristic time for the system to heat by one 
phonon starting from the ground state is T g r^huj m /^gk B T . Note that r" 1 does not necessarily 
reflect the actual collision rate between the sphere and gas molecules, R co \\~it Pvr 2 / k B T (it 



is possible for a single collision to be quite rare, Rcoii^Tg 1 , and to impart several phonons 
at once). We define a mechanical quality factor Q g = u) m /^ g due to the background gas, and 
a number of coherent oscillations Noil = uj m T g /2Ti expected before the energy increases by a 
single phonon. For a sphere of radius r = 50 nm, u) m j (2tt) = 1 MHz, and a room-temperature 
gas with P = KT 10 Torr, one finds 7 9 ~1(T 6 s" 1 ,Q 9 ~6xl0 12 , iV^J^lO 5 , indicating that the 
levitated sphere can be essentially de-coupled from its thermal environment. 

Photons scattered by the sphere out of the cavity lead to heating via momentum re- 



coil kicks. In analogy with atoms or ions trapped in the Lamb-Dicke regime [21| (when the 
particle is trapped on a scale Ax much smaller than A), the scattering induces transitions be- 
tween consecutive harmonic oscillator levels n— mil, with rates R n ^ n ±\ = 7 sc (n+l/2±l/2). 



Second-order perturbation theory 



30] yields 



7 SC = (2/5)(u r /uj m )R sc , (2) 

where u r = hk 2 /2pV is the recoil frequency and R sc = 487r 3 p^(^) 2 is the photon scatter- 
ing rate. A result identical to Eq. ([2]) holds for a weakly excited, trapped atom [31] . When 
photon scattering dominates the heating, the expected number of coherent oscillations is 

N (sc) = = 5 6+ 2 A 3 

osc 2tt 7sc 8tt 3 e - 1 V ' { } 

Note that ivisc scales inversely with the sphere volume (iV"isc^~40 for r = 50 nm, A = 1 /im, 
e^l), due to the fact that the scattered power and dipole force scale like pf nd and pi na \, 
respectively. Comparing with background gas collisions at P = 10~ 10 Torr and u m /(27t) = 
1 MHz, recoil heating dominates N osc for sphere sizes r^5 nm. Reaching the regime iV osc ;^>i 
implies that the sphere can coherently evolve for many oscillation periods after any cooling 
mechanisms are turned off, which makes this system a promising candidate for observing 
coherent quantum phenomena. 

Finally, we remark that R sc can be very large (i? sc ~10 15 s -1 for I = 1 W//im 2 and 
r = 50 nm) compared to atoms or ions, which enables direct imaging. The large scattering 
is due to the large intensities and the linear response of the sphere (it is not saturated like 
an atom or ion), as opposed to the system behaving as a lossy element in the cavity. The 
contribution to the cavity loss rate is k sc = 12tt 2 lj(V 2 /X 3 V c ) , where V c is the cavity 

mode volume, and is typically much smaller than the natural cavity linewidth k. We also 
emphasize that in the Lamb-Dicke regime, the scattering does not cause extra decoherence 



beyond that from recoil heating. This is in contrast to motional wavepackets of spatial 
extent Ax ~ A, where a single scattering event can destroy quantum coherence 32]. 



II. COOLING THE CENTER-OF-MASS MOTION TO THE GROUND STATE 

We now describe the optical cooling effect of the weaker, second cavity mode (denoted 
mode 2). For concreteness, we assume that the sphere is trapped near the anti-node x = of 
cavity mode Ei<x cos kix, and that the second mode has spatial profile E 2 cc cos (k 2 x — vr/4), 
such that the intensity gradient is maximized. The total Hamiltonian of the system is given 
in a rotating frame by 

H = —%5\d\d\ — h5 2 d\d 2 H — — (Si + a\) + \/2C ! '{d 2 + d 2 ) 

—hgi(cos 2k\X — l)a\ai — Tig 2 cos 2{k 2 x — 7c/4)a 2 a 2 + - — . (4) 

Here p, x are the momentum and position operators of the CM, di is the photon annihilation 
operator of cavity mode i, and Q, Qy/2£' are the driving amplitudes of modes 1 and 2, 
respectively. 5i is the detuning between the driving field and mode frequency when the 
sphere sits at x = 0. The opto-mechanical coupling strengths gi = ^—j^^i characterize 
the position-dependent frequency shifts due to the sphere (see Appendix), where V c ,i, uii are 
the mode volume and resonance frequency of mode %. To simplify notation, we assume that 
modes 1,2 have similar properties, u;i~u;2 = oo, etc. In addition to the evolution described 
by H, the system also exhibits cavity losses and the mechanical noise described previously. 

Expanding the opto-mechanical coupling term of mode 2 around x = 0, hg cos 2{kx — 
Tr/4)a\a 2 ~ 2hgkxa\a 2 , one finds a linear coupling in the sphere position, analogous to the 
effect of radiation pressure on a moving mirror of a Fabry- Perot cavity [lO] . Physically, the 
motion induces changes in the detuning and thus the intra-cavity field amplitude, while the 
lag in the cavity response enables the field to do work (cooling) on the sphere. To calculate 
the cooling rate, following the techniques of Ref. 1C| we first apply shifts to the operators, 
di — > di + «j, x — > x + x , where «j and x ~ (/k (( w k 2 ('/(k 2 + 45|)) are mean values 
of the quantum fields. Here we have defined 2( = \a 2 /ai\ 2 as the ratio of intra-cavity 
intensities of modes 1 and 2, and assumed that mode 1 is driven on resonance (5i = 0). To 
lowest order in £, field mode 1(2) is purely responsible for trapping (cooling). Subsequently 
tracing out the cavity degrees of freedom yields equations for the mechanical motion alone. 



In particular, to lowest order in ( and for 5 2 < 0, the cooling laser provides a net cooling 
rate T = R opt ,_ - R opt , + = Kn 2 m [((5 2 + u m ) 2 + (k/2) 2 )- 1 - ((5 2 - u m ) 2 + (k/2) 2 )' 1 } (see 
Appendix), where R op t^ denote the anti-Stokes (cooling) and Stokes (heating) scattering 
rates (see Fig. [Ib). Here Q m = 2gkx m \al\^/2( is the effective opto-mechanical driving 
amplitude (see Fig. [lb) and x m = ^Jh/2muj m . Validity of these perturbative results requires 
Q m < K,u m and £<1. 

In the realistic limit that background gas collisions are negligible, the steady-state phonon 



number is (n/)~ny+7 sc /r, where fif = R opt)+ /T is the fundamental limit of laser cooling [10 ]. 
It is minimized when 8 2 = — (1/2) ^/k 2 + 4^. In particular, when sideband resolution is 
achieved (u m > k), nj min ^(K/4u; m ) 2 <^l, indicating that ground-state cooling is possible 
provided other heating mechanisms are made sufficiently small. Considering the limit uj rri ^>n 
and taking the maximum cooling rate r~K consistent with the perturbative calculations, 
using Eq. ([3]) one can then re- write (n/) as 

<"/>«l^ + ^- ( 5 ) 

The last term on the right corresponds to photon recoil heating and (ft = (47r 2 /5)(V7A 3 )^ 
is a dimensionless parameter characterizing the sphere volume. Eq. ([5]) is minimized for 
K,/uj m = 20 1 / 3 , in which case (%} j = 3(/> 2 / 3 /4oc(r/A) 2 ^Cl. Thus, one sees that ground-state 
cooling is readily attainable (provided that can be simultaneously satisfied). Physically, 
the optimum value of n/uj m balances good sideband resolution and excessive heating when 
large intensities are used to increase u m . 

To illustrate these results, we consider a sphere of radius r = 50 nm and io m j (2tt) = 
0.5 MHz levitated inside a cavity of length L = 1 cm and mode waist w = 25 /im (V c = 
(tt/4)Lw 2 ). In Fig. [2^ we plot the minimum obtainable (rif) (black curve) as a function of 
cavity finesse T = itc/2kL, assuming negligible gas collisions and subject to the constraints 
(,n m / K,Q m /uj m < 1/2 and optimized over detuning 5 2 - For low cavity finesse the cooling 
is essentially limited by sideband resolution (n/ min , red curve) and the ground state regime 
(nf)r^l can be reached with a finesse of JF~3600. A minimum of (n/)~0.02 is reached at 
a finesse of JF~50000 (with corresponding cooling rate r~10 6 s -1 ). This corresponds to a 
final temperature of Tj~6 /zK, or a remarkable compression factor of T/Tj~5xl0 7 relative 
to room temperature T. 



7 



III. MOTIONAL ENTANGLEMENT AND SQUEEZED LIGHT GENERATION 
USING QUANTUM STATE TRANSFER 



A number of related schemes have been proposed for quantum state transfer between 



light and the motion of atoms 
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13] or nano-mechanical systems \l4 . Il5| . In our system, 
the small mechanical noise and ease of achieving good sideband resolution in principle allow 
state transfer to be realized with almost perfect efficiency This might enable light with 
non-classical properties to be mapped onto mechanical motion, and as an example, we show 
that this can be used to generate EPR correlations between two spatially separate spheres. 
Moreover, a complementary process can be realized, where a non-trivial mechanical state (a 
squeezed state) is prepared through coherent manipulation and subsequently transferred to 
light leaving the cavity. The latter case illustrates how opto-mechanics can yield a novel 
nonlinear optical system. 

First we give a simplified picture of quantum state transfer using a one-sided, ideal 



cavity (where all losses are via transmission through one cavity mirror) |33|]. Specifically, 
we consider the Heisenberg equations of motion in a rotating frame for the cavity cooling 
mode and the motion (after applying the shifts described in the previous section), when the 
cooling mode is driven resonantly on the red motional sideband (5 2 = —co m ), 

—a 2 = --a 2 - itt m (b + tfe 2iuJmt ) + yfnd^ 
jb = (i/h) [H e , b] - in m (a 2 + 4e 2 ^ *) + iF(t)e^. (6) 

The Hamiltonian H e describes any external forces or couplings applied to the sphere beyond 
those in Eq. (jl]), b is the annihilation operator corresponding to a harmonic oscillator of mass 
m and frequency u> m , and d 2 ^ n is the cavity input operator associated with losses. F(t) is 
the (Hermitian) noise due to photon recoil, which has correlations (F(t)F(t')) = <j)u m 5(t—t'), 
and we assume all other noise is negligible. Since the cavity trapping mode (d\) effectively 
provides a harmonic potential and can otherwise be ignored, for simplicity we will omit the 
subscript 2 as we refer to the cooling mode in future discussions. Temporarily assuming that 
the non-secular terms {e 2lu)mt ) can be ignored and that the mechanical motion evolves slowly 
on time scales compared to l//t, one can adiabatically eliminate the cavity mode to yield 
d^ - 2i{Q m /n)b + (2/VK)a in , and db/dt^(i/h) [H e , b] - {T/2)b - i^fd in + iF{t)e iuJmt , where 
r=4f2^//€ is the cavity- induced cooling rate in the weak-driving limit (Q m ^ k). The cavity 
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output is related to the input and intra-cavity fields through a out = ^/na — a in [33|, or a ut~ — 
i\/Tb + a; n , which states that the mechanical motion is mapped onto the outgoing light. 
Physically, the cooling process converts phononic excitations into photonic excitations that 
leak out of the cavity. Generally, two mechanisms will degrade state transfer. First, F adds 
extra noise to the ideal state that one is attempting to transfer, with a strength characterized 
by the small parameter 0. Second, the non-secular terms contribute to Stokes scattering, 
destroying the perfect phonon-photon correspondence, with a strength that is expected to 
be proportional to (n/oj m ) 2 . Given that 0, (K,/uj m ) 2 can be made small, nearly perfect state 
transfer is possible in principle. We illustrate this with two examples, entanglement transfer 
and squeezed light generation. 



A. Entanglement transfer 

Here we describe how EPR correlations shared between two modes of light [17[ can be 
mapped to the motion of two spheres trapped in spatially separate cavities. Specifically, we 
define quadrature operators for the input light for each of the two systems (denoted A, B), 
given by X+\ n = (a^ + a[^), X^\ n = (aj^ —a[^)/i for j = A,B. A similar set of operators 
Xf$P m , X± out can be defined for the motion and output light, by replacing a^— >W\ a^ t , 
respectively. Of particular interest is the case where the two input fields exhibit broadband 
EPR correlations between them, 

<(4iM + *&(<")) 2 >/2 = ((X™(u>) - X^)f)/2 = e~ 2R < 1. (7) 
When the variances satisfy e~ 2R < 1, the two modes exhibit correlations below vacuum level 



and are entangled 34| (for concreteness, we assume the other combinations of quadratures 



satisfy ((I^(w)Tlf i ! 1 (w)) 2 >/2 = e 2R ). Such EPR correlations have been observed with 
light and in the internal degrees of freedom of atomic ensembles 35|, but have yet to be 
demonstrated using mechanical systems. 

To proceed, we solve Eq. (J6j) in the Fourier domain (including the non-secular terms) for 
each of the systems for the correlations given in Eq. (J7J) and H e = 0. Generally, the non- 
secular terms yield an infinite set of algebraic equations (coupling frequencies u m + 2nuj m 
for integer n), which given uu m ^>n m ,K can be truncated to good approximation at n > 1. 
For simplicity of analysis, we assume the two systems have identical properties, and that 
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the cooling rate r ~ k. However, we expect our results should qualitatively hold provided 
that only T, u m of the two systems are properly tuned, which can be easily accomplished by 
adjusting the trapping and cooling beam intensities. One can then show that state transfer 
yields the following joint variances in the motion (see Appendix), 

A EPR ee ((li!»W T lSW) 2 )/2 = e~ 2R + T^(3e 2 * + 2sinh2i?) + (8) 

As expected, Stokes scattering and recoil heating contribute to the variance by amounts 
(K,/oj m ) 2 and <pu m /K, respectively. This can be minimized with respect to n/uj m , yielding 
A E pR,min = e _2fi + 3(0/2) 2 / 3 (3e 2i? + 2 sinh 2R) 1 / 3 . To illustrate these results we plot A EPR>min 
in Fig. [2b as a function of e~ 2R , taking the same parameters as in Fig. [2^. For the moderate 
values of e~ 2R typically obtained in experiments |17|. EPR correlations in the motion can 
be achieved with reasonable cavity finesse F < 10 5 . 

B. Squeezed light generation 

First we describe a technique to create a mechanical squeezed state, and then derive 
the properties of the outgoing light upon quantum state transfer. Mechanical squeezing is 
accomplished by adding a sinusoidally-varying component to the intensity of the trapping 



beam, which yields the Hamiltonian of a parametric amplifier 36( , H e = emUJ^x 2 sin 2uu m t. 
Here e m is a small parameter characterizing the strength of the modulation of the trap fre- 
quency. As one approaches the threshold for parametric oscillation (e m co> m — ►r), the variance 
in one quadrature of motion is reduced by up to a factor of 2 36]. 

We now investigate the properties of the outgoing light over a narrow frequency range 
near the cavity resonance, specifically considering X± )0nt (uj = 0). We apply similar methods 
as above to solve Eq. ([6]) in the Fourier domain. Taking the limit as one approaches threshold 
and the variance in the output light is given by (see Appendix) 



AX + , out (u, = 0) = — + --. (9) 



Again, an optimum value of k/ uj m (xcf) 1 ^ maximizes the squeezing, with 
(AA 2 out ) min ^2.O40 2 / 3 (note that AX 2 out = 1 for vacuum). A plot of (AX 2 out ) min 
as a function of sphere size is shown in Fig. [2b. For r = 10 nm size spheres, one finds that 
over 25 dB of squeezing relative to vacuum can be obtained using an ideal cavity (note 
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for good vacuum conditions, background gas collisions are negligible down to r~5 nm). In 
practice, a cavity has additional scattering and absorption losses that limit the squeezing. 
Taking an ultra-high finesse cavity with ~1 ppm losses per round trip js^] and a set of 
reasonable cavity dimensions, we show in the Appendix that light squeezed by up to ~15 dB 
can be extracted. 

In principle, similar techniques also apply to trapped atoms or ions. However, one benefits 
from the relatively large mass m of the sphere. Specifically, approaching threshold, one 
quadrature of motion becomes infinitely unsqueezed, producing a large position uncertainty 



Ax |36|]. At the same time, faithful quantum state transfer requires a linear opto-mechanical 
coupling, which translates into a requirement that the Lamb-Dicke parameter i]=kAx oc 
m~ l l 2 <^l be small. In the Appendix, we show that ^-ClO -2 can be satisfied with a sphere 
even in the regime of ~25 dB squeezing levels. To compare, a typical atom trapped with a 
frequency of c<j m /(27r)~l MHz in its ground state already yields r]~0.05. 

IV. OUTLOOK 

An optically levitated opto-mechanical system can have remarkably long coherence times, 
which potentially enables quantum phenomena such as entanglement to be observed even 
in room-temperature environments. Combining previously demonstrated techniques to con- 



trollably grow small particles [38( and load and manipulate them in vacuum p, |39( should 
put this system within experimental reach. Extending the ideas presented here should open 
up several other interesting possibilities. For example, beyond the dipole-type objects con- 
sidered here, one could engineer the shapes of the levitated objects to yield even larger 
mechanical frequencies and coherence times, and controllably study the decoherence of a 



large system [40|. Also, several spheres or more complex nano-mechanical systems with in- 
ternal modes could be levitated and coupled together, for the purpose of entangling multiple 
degrees of freedom. Separately, one could take advantage of the fact that the potential for 
the CM is purely optical to engineer non-trivial dynamics, such as nonlinear motion. It 
would also be interesting to develop analogous levitation techniques using nano- and micro- 



photonic cavities 4jJ, |42|, combining their remarkable optical properties with the excellent 



mechanical characteristics of a trapped particle. Finally, by levitatin g ch arged or magnetic 



systems, one could potentially realize systems analogous to ion traps 43( or facilitate novel 
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quantum hybrid architectures 
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APPENDIX A: ELECTRODYNAMIC CALCULATION OF FORCES ON DI- 
ELECTRIC SPHERE 

Here we describe in detail the point-dipole approximation for a sub-wavelength dielectric 
sphere, and compare the results of this approximation with exact numerical electrodynamic 
calculations of the optical forces. For concreteness, we consider a dielectric sphere of (pos- 
sibly complex) permittivity e and radius r, interacting with an incident standing electro- 
magnetic wave with electric and magnetic field components E; n = fi^cos k(z — Zq) cos ut 
and B in = y(E /c) sin k(z — z ) sin ut (k = u/c). Here we assume that the sphere is in 
free space rather than in a cavity, which allows one to unambiguously calculate the optical 
forces acting on the sphere independent of its motion (as opposed to the cavity case where 
the motion of the sphere generally shifts the cavity resonance, causing the intra-cavity field 
Eo(t) to depend on the history of motion). In the special case where the sphere is local- 
ized near one of the nodes or anti-nodes, the free-space and cavity cases yield the same 
results (e.g., for the mechanical trap frequency uj rn ) as the cavity resonance and intra-cavity 
field to lowest order become insensitive to small displacements of the sphere. The electrody- 
namic problem of plane- wave scattering off of a sphere is exactly solvable, as the vector wave 
equation V 2 E(r) + /c 2 e(r)E(r) = (with similar equation for B) admits solutions through 
separation of variables [451 ] . Note that one can define natural length scales k\y/e\r, kr for 
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the electrodynamic response inside and outside the sphere. Of particular interest is the case 
when fc|\/i|r<Cl is a small parameter (we assume that > 1 for this discussion, which 
is typically the case). One can then formally solve the wave equations using perturbation 
theory, with the lowest order equation given by V 2 E(r) = along with appropriate bound- 
ary conditions at the surface of the sphere. Physically, this approximation states that the 
magnetic field is not important in the near-field, such that the lowest-order response of the 
sphere can be obtained by solving an electrostatic equation. Taking an optical wavelength 
of A = 2n/k = 1 /iin and e = 2, for instance, the electrostatic solution should be valid for 
r<110 nm. In this regime, the polarizablity of the sphere is of the simple form given by 



electrostatic theory, ai n d = 3eo^fr| |45| . as is used in the main text. The optical poten- 
tial experienced by the sphere is predicted to be U opt = — (l/4)(Re ai n d)-Eo cos 2 k[z — z ). 
For spheres larger than r>l//c|\/e|, the forces predicted by the electrostatic theory will be 
substantially larger than the actual forces, as phase variations of the field within the sphere 
become important. 

To compare the electrostatic approximation with actual results, we first solve the elec- 
trodynamic scattering problem exactly. The exact force F z along z can then be obtained by 
integrating the Maxwell stress tensor over the sphere surface S, 

F z = e j> da ^ T zj^j, ( Al ) 

" j=x,y,z 

where hj is the outgoing normal vector to the sphere surface. In Fig. [3] we compare the 
approximate and exact forces for various values of r, taking e = 2. It can be seen that the 
two methods agree closely for k^/er^l. For spheres where ky/er>l, the forces predicted 
from electrostatic theory can be much larger than the actual forces, and even different in 
sign. 



APPENDIX B: ABSORPTION LOSSES OF TRAPPED SPHERE 

In this section, we consider the effect that a small imaginary component of the permit- 
tivity e has on a trapped sphere. In the limit that the sphere has a radius much smaller 
than the optical wavelength, the sphere behaves as a point-like dipole with polarizability 
«ind = 3eoV (fr|j- For small Im e, the polarizability acquires a small imaginary component 



13 



that leads to a non-zero absorption cross-section, with a corresponding absorbed power 

P abs = 127r-Vlm (Bl) 
A e+2 K ' 

Here Jo is the trapping beam intensity, V is the sphere volume, and A is the optical wave- 
length. The absorbed power causes a rise in the internal temperature T int of the sphere, 
which is balanced out by thermalization with a background gas and blackbody radiation. 

We first quantify the effect of the background gas (which is negligible in the regime of 
particular interest where the sphere is trapped under good vacuum conditions). There are 
two limiting regimes to the background gas interactions, where the sphere radius is either 
much smaller or larger than the molecular mean free path A m f p . At a relatively large pressure 
of P = 1 Torr and room temperature, A m f p ~100 f/m and thus our case of interest is always 
r^cA m fp. Here, gas molecules independently collide and partially thermalize with the sphere. 



This leads to a cooling rate [46( 



~dt = - a »y^ {lTr )Pv -^r {-t ~ l ) ' (B2) 

where P, v TmB , T are the background gas pressure, root-mean-square speed, and temperature, 
respectively, and 7 s h is the gas specific heat ratio (7 s h = 7/5 for an ideal diatomic gas). a g 
is a phenomenological energy accommodation factor (0 < a g <l), which characterizes the 
degree to which a gas molecule thermalizes with the sphere upon a single collision. 

Under good vacuum conditions, blackbody radiation dissipates the majority of the power 
absorbed by the sphere. For the sub-micron spheres we are considering, the radius is much 
smaller than the absorption length at typical blackbody radiation wavelengths, and thus 
the usual formulas for blackbody radiated power do not apply. Instead, the sphere again 
behaves as a point-like dipole at these wavelengths, e.g., the radiated power scales like 
volume (as opposed to surface area in the case of a large object). The internal heating rate 
due to blackbody radiation is given by dE/dt = Xlk(^ c ^)-^ abs > k ' wnere the sum is over all 
blackbody radiation modes (and polarizations), k is the wavevector of each mode, and -R a bs,k 
is the absorption rate of each mode. It is given by 

/W = 3ck(V/V q )n k lm ~ ^ , (B3) 

where n k = {e %ck / kBT -I)' 1 is the occupation number of each mode and V q is the quantization 
volume. Assuming that the sphere has a relatively constant and temperature-independent 
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permittivity e(u)^iehb across the blackbody radiation spectrum, it is straightforward to show 
that the sphere absorbs blackbody radiation at a rate 

dE - 72C < 5 > V Am (k B T)\ (B4) 



dt 7T 2 c 3 h 4 V e bb + 2 / 

where T is the background temperature and £(5)~1.04 is the Riemann zeta function. Simi- 
larly, the sphere radiates blackbody energy cit £t rstte given by the negative of Eq. flB4j) . with 
the substitution T— >T int - 

To illustrate these results, in Figs. H^-c we plot the internal equilibrium temperature T int 
of the sphere as a function of background gas pressure and trapping intensity Iq. Here we 
have taken into account the effects of optical absorption (Im e), thermalization with the 
background gas, and blackbody radiation. The values of Im e in Figs. H^i,b,c correspond to 
bulk optical absorption rates of 10, 100, 1000 dB/km, respectively, while the real part of the 
permittivity is chosen to be Re e = 2. We have taken the other parameters to be r = 50 nm, 
a g = 0.25, Im = 0.1 (roughly corresponding to the averaged value of fused silica around 
blackbody wavelengths [47J), and a volumetric heat capacity of the sphere of c = 2 J/m 3 -K. 
Note that at sufficiently low pressures, the temperature becomes pressure-independent as 
only blackbody radiation significantly contributes to energy dissipation (as indicated by 
the vertical contours in the figure). Furthermore, in this regime the final temperature is 
independent of the sphere size (provided that r<^cA), since both the optical absorption and 
blackbody radiation scale linearly with volume. For losses of ~10 dB/km, one finds that 
over 10 W//im 2 of power can be sustained without exceeding the melting point of a typical 
material. 



APPENDIX C: DERIVATION OF OPTO-MECHANICAL COUPLING 
STRENGTH 

Generally, introducing a dielectric material into an optical cavity causes the bare resonant 
frequency a; of a cavity mode to shift by an amount 5u, which in perturbation theory is 
given by [48( 

5u 1 / rf 3 r 5P(t) ■ E(r) 

u 2 / rf 3 r e E 2 (r) ' 1 ' 

Here E(r) is the bare cavity mode profile and 5P(r) is the variation in permittivity introduced 
by the dielectric object. Considering the case where the dielectric object is a sub- wavelength 
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sphere, its dielectric response is well-approximated by a point dipole, P(r')&ai n dE(r)5(r — 
r'), where r is the center-of-mass (CM) position of the sphere. Taking a mode profile 
.Eoccos (kx — 0), one readily finds (up to a constant shift) that 

3 €. 1 

8u) = - — — cos (2kx - 2(j))u. (C2) 

The interaction Hamiltonian between this optical mode and the mechanical motion is sub- 
sequently given by H om = HScuaJa, and as in the main text, one can define a characteristic 
opto-mechanical coupling strength g = f^rf^^- 

APPENDIX D: OPTICAL SELF-COOLING EQUATIONS 

Here we derive in detail the cooling rate equations for the CM motion of the sphere, whose 
results are summarized in the main text. We begin with the Hamiltonian given by Eq. (4) 
in the main text. The corresponding Heisenberg equations of motion, including dissipation, 
are 

^j-ai = (iSi - re/2)ai - y + y/Hai^, 
^-a 2 = (i(5 2 + 2gkz) - n/2)a 2 - ^^2^ + V*>a 2 ,in, 
—p = — 4hgk 2 a\aiz + 2hgka\a 2 — 7^/2 + F p (t), 

(JjL 

U = l - (Dl) 
dt m v ' 

Here are input-field operators associated with the cavity mode losses k, 7 is the damping 
rate of the motion, and F p is the noise force acting on the sphere. In the above equations, 
we have expanded the position-dependent opto-mechanical coupling terms gi cos 2{kiX — (pi) 
to first order in the displacement x, and for simplicity have assumed that the two cavity 
modes have similar properties (gi~g 2 = g, etc.). We now apply shifts to all of the operators, 
ai^cii + ai, x— >x + xo, where the constants xq and ctj are chosen to cancel out all of the 
constant terms in the equations of motion. This yields 

ai = --, (D2) 

K 



a2 = "T(«/2)-^' (D3) 
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where 8' 2 = 62 + 2gkxa is the detuning relative to the new resonance frequency of the cavity 
when the sphere sits at x = xq rather than x — 0. Physically, x = xq corresponds to 
the minimum of the total optical potential formed by the two driven cavity modes. We 
define the ratio of the cavity mode intensities to be 2( = |a2/«i| 2 , which is equivalent to 
( = ('k 2 /(k, 2 + 4<5 2 2 ). In terms of (, the shifted equilibrium position is given by kxo = (. 
Clearly then the expansion in x of the opto-mechanical coupling terms requires that ( be 
small. For simplicity, the prime symbol in 5' 2 will be implicitly understood, and we also take 
5\ = in the following discussions. Following the shifts to the operators and x and then 
linearizing the equations of motion, one finds 

— ai = — Aigk 2 XoaiX — (/t/2)di + \A^i,in> 
^-a 2 = {iS 2 - n/2)a 2 + 2iga 2 kx + a/k^,^, 

— p = —Ahgk 2 \ai\ 2 x + 2Kgk ^02^2 + a 2^2 — 2kxo(aia\ + a : ai)^ — 7p/2 + F p (t), 

—x = p/m. (D4) 
at 



Note that cavity mode 1 provides a linear restoring force dp/dt^ — Ahgk \a±\ x = —mtu^x, 
and it is straightforward to show that this relation leads to the expression for the harmonic 
oscillator frequency u m given in Eq. (1) of the main text. Furthermore, note that the sphere 
is opto-mechanically coupled to mode 1 with an amplitude 4gk 2 x ai(x(, and to mode 2 with 
an amplitude 2gka2(x^f(. Thus, to lowest order in (, modes 1 and 2 are purely responsible for 
optical trapping and cooling, respectively. Treating mode 1 simply as an external harmonic 
potential for the sphere, the opto-mechanical system comprised of the CM motion of the 
sphere and cavity mode 2 is completely equivalent to the system described in Ref. [10 ] . In 
particular, the optical self-heating and self-cooling rates R± given in the main text follow 
immediately. For convenience, we also re-define the phases of the operators to make the 
opto-mechanical driving amplitude Q m = 2hgka>2 = 2hgkai^/2( real. 



APPENDIX E: NOISE FORCES ACTING ON TRAPPED SPHERE 

In the main text, we have derived the motional heating rates of the sphere due to back- 
ground gas collisions and photon recoil kicks, which under realistic conditions are the dom- 
inant heating mechanisms. Here, we derive the heating rates for a number of other less 
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important processes. 



1. Photon shot noise 

Photon shot noise inside the cavity leads to heating via fluctuations in the mechanical 
oscillator frequency u m . We write the varying mechanical frequency in the form 



<(t) = <o \l + J > ( E1 ) 

where u>m t o, N are the mean frequency and mean photon number in the trapping mode of 
the cavity, and 5N is the number fluctuation of this mode. Following the techniques of 
Ref. {49J , the shot noise leads to parametric transitions (where the phonon number n^n±2 
jumps in pairs) at a rate R proportional to the power spectral density of the fluctuations at 
frequency 2u mfi , 

R n -, n+2 = —^-S{2u m , ){n + 2 n + 1 , (E2) 
lo 

Rn^n-2 = ^ ™'° S(2uj m , )n(n- 1). (E3) 
lo 

Here the power spectral density is defined by 

S(u) = 5 / dt cos uot (5N(t)6N(0)), (E4) 

which is evaluated to be S(u) = ^2^.4^2 for a cavity of linewidth k driven on resonance. 
Assuming that the sphere initially is in the ground state, the number of oscillations before 
a quantum jump due to shot noise is 

Ar(sn) U rnfi 6 + 2 V c p L0 m o , 2 , 2 \ /T?r-\ 

NL = 2^7 2 = 7=i*rt&— {K + 16 "- o) - (E5) 

Here, k = 2tt/X is the wavevector of the trapping beam and V c is the cavity mode volume. As 
an example, we consider a cavity of length L = 1 cm and waist w = 25 /xm (V c = (it/4)Lw 2 ), 
a high-index dielectric sphere (^r|~l), density p = 2 g/cm 3 , A = 1 pm, and trapping 
frequency co> m /(27r) = 0.5 MHz. A^osc^ as a function of cavity finesse F (F = ttc/2kL) is 
plotted in Fig. [51 It can be seen that the number of allowed oscillations is at least of order 
A^c^~10 10 , which is much larger than the limit due to photon recoil. Physically, the low 
heating rates are attributable to the large intra-cavity intensities used to achieve ~MHz 

— 1/2 

mechanical oscillation frequencies, which suppresses the fractional noise 6N/N ocN 

18 



2. Blackbody radiation 



As in the case of scattering of laser light, the absorption and emission of blackbody radi- 
ation by the sphere also lead to recoil heating. The absorption rate of blackbody radiation 
of mode k is given in Eq. (1B3[) (with each absorption event providing a momentum kick hk x 
along the trapping axis), and again we assume that e(u;)~e(,b is approximately flat across 
the blackbody radiation spectrum. Summing over all modes, the characteristic jump rate 
due to absorption of blackbody radiation is then given by (cf. Eq. (2) in main text) 

2vr 4 {k B Tf e bb -l ( 

7bb = -^-^75 Im (E6) 

bS a>n puj m e bb + 2 

The jump rate between harmonic oscillator levels is R n ^ n ±i = 7 bb (n + 1/2=1=1/2). An analo- 
gous expression holds for heating via the emission of blackbody radiation, with the replace- 
ment T— >Tinf Note that 7^ is size-independent for small spheres, as both the absorption 
rate and mass scale linearly with V. Taking as an example a system with u; TO /(27r)~l MHz, 
p = 2 g/cm 3 , Im ^ bb ~ 2 = 0.1, and T~T int ~300 K, we find that the number of oscillations 
before a quantum jump (due to either absorption or emission) is ivjgc ~10 n . 



3. Anisotropy of sphere 

The general problem of the rotational motion of an arbitrary dielectric object inside an 
optical cavity is quite challenging to solve. Generally, the polarizability becomes a 
function of its orientation, and changes in its orientation lead to changes in the optical 
trapping potential and the intra-cavity intensity. Here we consider a simplified version of 
the problem, where the rotational motion is limited to one axis, and the anisotropy or 
deformation of the sphere is of spheroid- type. As in the case of the sphere, the latter 
assumption admits analytical solutions for the polarizability tensor of the object [50]- In 
particular, we assume that the dielectric is a prolate nanospheroid whose size is much smaller 
than the optical wavelength, with semi-major axis a and semi-minor axis b, and that the 
ratio a/bttl [i.e., the deviation from an ideal sphere is small). Then the polarizability of 
the spheroid is given by 

«ind^« indi0 (l±^^ [(«/&) 4/3 " 1]) (E7) 
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with amd,o~3eo^f+2. Here the ± symbols denote when the major and minor axes are aligned 
along the field polarization axis, respectively. From Eqs. fIClj) and flE7j) . it is straightforward 



to find the shift in the cavity frequency taking into account the rotational degree of freedom, 

5uj = 5u + 5uj0 cos 29, (E8) 

where 5u is the shift associated with the CM position alone (as given by Eq. ( 1C2I) ). and 

27 V f e — 1 \ 2 

J [(a/6) 4 / 3 -l]^cos(2fcr-20). (E9) 

Here we have defined 9 as the angle of rotation of the spheroid. 

We are now interested in deriving the effect of the rotational motion on the CM motion. 
In analogy with Eq. ( ID If) . the coupled equations of motion between the rotation and the 
trapping mode are 



— (iduje cos 29 + — j a± + 



dt 

= 2h5uj e \ai\ 2 sin 29 - >y e p e + F e (t), 

dt 

d9 pa s 

* = V (E10) 

where pg is the angular momentum associated with 9, Ig is the moment of inertia, and 
70, Fq are the damping coefficient and noise force acting on the rotational motion. Since 
the rotational energy is of order ~fceT, it suffices to consider the classical equations given 
here. The damping term is effected through the background gas, as each collision between 
the spheroid and a gas molecule partly exchanges angular momentum between the two 



systems. The damping coefficient is found to be 70 = 5 a/3/ (2ir)agP/ (v rms rp) [5lJ, where 
r^a^b. ag is a phenomenological accommodation coefficient describing the efficiency of 
angular momentum transfer. The noise force has correlations (F(t)F(t')) = 2D5(t — t'), 
where D = joksT '/Ig. Note that 70 is a very small quantity under good vacuum conditions. 

The full nonlinear coupled equations of Eq. (1E10I) are difficult to treat in a general set- 
ting. However, given the typical smallness of the parameters 5uog/n and hSug\ai\ 2 / '{ksT) 
for nearly spherical particles, to lowest order we can ignore the optical coupling to the ro- 
tational motion, and the dominant effect of the sphere anisotropy is trap heating through 
fluctuations in the polarizability ct in d rather than intra-cavity intensity fluctuations. This 
leads to fluctuations in the trap frequency given by 

Su m (t) = €gUJ mfi cos 29(t), (Ell) 
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where eg = -§)j^(( a /b)^ 3 — 1). In analogy with the discussion in Sec. IE 1[ these fluctuations 
lead to parametric heating, with a jump rate out of the ground state given by 

poo 

Ro^2= dt cos 2u mfl t (8u m (0)8u m (t)}. (E12) 

J 

Denoting 59 (t) = 9(t) — 9(0), the above equation can be re- written in the form 



OO 



Ro^2 = \j dt cos 2cu mfl t (e e a4 ) 2 (cos 289{t)). (E13) 

Making a Gaussian approximation (e 2 * <5e ^- ) )~ exp(— (59 2 (t)) /2), and taking the limit of small 
70, one finally finds 

^ = <8 M5y exp V m)- (E14) 



Here u r = d6/dt is the angular velocity of the spheroid (typical values of a/ (uj 2 ) are in 
the MHz range for sub- wavelength particles). Note that the above function is peaked at 
uj m fi = \/ (on 2 ), i.e., the parametric heating is most pronounced when the rotational fre- 
quency is comparable to the CM oscillation frequency. At this maximum, i? ^ 2 /^m,o~0.2eg. 
Furthermore, for this worst-case scenario, Ro^/^mfl can be suppressed to the ~1CT 5 level 
with an anisotropy of a/6~1.03. 

APPENDIX F: ANALYSIS OF ENTANGLEMENT TRANSFER 

Here we provide a detailed analysis of entanglement transfer between two modes of light 
and two spatially separate spheres, leading to Eq. (8) in the main text. The EPR correlations 
between the two light modes given by Eq. (7) in the main text, 

((4lM + ^?») 2 >/2 = ((*<&M - x[%(u)) 2 }/2 = e- 2R < 1, (Fl) 



are of the form created by a non-degenerate optical parametic amplifier (NOPA) 52] , which 
we describe below. 

The Hamiltonian corresponding to a NOPA with cavity modes A, B is given by 

H = ih(J3/2)(cWcW - c^c™), (F2) 



where is the annihilation operator of mode j. Taking an ideal, one-sided cavity [33j, and 
assuming that the modes have identical linewidths k c , the Heisenberg equations of motion 
for each mode read 



dt ^ = -^-^ + ^l>- (F3) 
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Here c in is the cavity input field for mode j, and f = A,B for j = B,A. The output 
field is related to the intra-cavity and input fields by c^ t = y/n^c^ — c^. Writing c^\t) = 

Eq. (1F3|) can be exactly solved in the Fourier domain for 
Specifically, defining quadrature operators X.^} = + and = (c^ — c^)/i (with 
analogous definitions for the quadrature operators of the input and output fields), one can 
show that 

Over a bandwidth Au;<Cfv c that is much smaller than the cavity linewidth, one can ignore 
the uj dependence in the equation above, yielding 

X { &M ± *£LM = e~ R (X^l(u) ± xW{uj) , (F5) 

where e~ R = for f3 < k c . Physically, for non-zero j3, the joint variance of these 

quadratures in the output fields can display reduced fluctuations relative to the input fields. 
It can also be shown that the other combinations of the quadratures (for Aco><^Ck c ) satisfy 

*J£» T XjgLiu) = e R (xW(u) =F xgiH) , (F6) 

such that their joint variances become enhanced. For this discussion, the input fields to the 
NOPA are assumed to be vacuum states. 

We now consider the quantum state transfer process for two spheres trapped in spatially 
separate cavities, where the two output fields generated by NOPA are fed as input fields into 
each of the opto-mechanical systems. The equations of motion for the two opto-mechanical 
systems (denoted A, B) are given by Eq. (6) in the main text, with the replacement = 
c~jrt. As in the main text, for simplicity we suppress the subscript "2" in the field operators 
denoting the trapping mode, since we are only interested in this mode from this point on. To 
solve these equations, we again work in the Fourier domain. Without the fast-rotating terms 
e 2luJmt , one could achieve ideal state transfer between the mechanical motion and light, as 
discussed in the main text. When the fast-rotating terms e 2luJmt are included in the analysis, 
the frequency components u,u + 2nu m (integer n) of the operators are coupled together 
in an infinite set of algebraic equations. To make the problem tractable, we truncate this 
infinite set by ignoring the components cft\u) + 2nuj m ),W\uj + 2ntu m ) where |n|>2 (e.g., 
we assume a^(u±4w m ) = 0). This truncation essentially amounts to the assumption that 
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uj m is large compared to the other frequency scales in the problem. We then solve the 
coupled set of equations for a^(uj),U j \uj) in terms of F^(uj) and a^(iv) (or c[^(u)), 
which allows us to obtain any correlation functions for the cavity field or mechanical motion 
in terms of those of the noise and input fields. The noise forces are assumed to be 
dominated by photon recoil heating and are independent for the systems A, B, such that 
their correlations take the form = and (F^(u)F^(u')) = <pu m 5(uj + u')5 jr , 

where = (47r 2 /5)(l / /A 3 )^| (see main text). We are specifically interested in the quantity 

A EPR ee ((xi%(t) T xi B l(t)y}/2 (FT) 

characterizing the joint variance in the motion of the two spheres. The solution is generally 
quite complicated, but can be expanded to lowest order in the small parameter n/u> m (it 
is reasonably assumed that sideband resolution can be achieved, so that K/u m <gX). After 
performing this procedure, and also ignoring any fast-rotating terms (e ±2luJmt ) in the final 
expression for A EPR , one arrives at the solution given by Eq. (8) in the main text. 



APPENDIX G: ANALYSIS OF SQUEEZED LIGHT GENERATION 

Here we derive the squeezing amplitude given in Eq. (9) of the main text. In the main 
text, it was argued that the trapping mode of the cavity can be effectively considered as a 
mechanical potential in the limit of small (. We consider the situation where the trapping 
beam intensity is varied to produce a sinusoidal component in the mechanical spring constant 
at frequency 2u m , with an effective Hamiltonian for the motion given by 
p 2 1 

H m = — + -muj 2 n x 2 (l + 2e m sin 2u m t) (Gl) 
= hu m ftb - ;M(£2 e W - StV 2 *""*) + 2 {hpftb sin 2cu m t} . (G2) 

In the last line, we have re- written x = \J^^~ (b + ft) and p = hm £ m (ft — b) in terms 
of the harmonic oscillator annihilation operator b and also defined (5 = euj m /2 (unrelated 
to the P term defined in the previous section for a NOPA). The term in braces is a fast- 
varying contribution to the Hamiltonian, in addition to the "ideal" squeezing Hamiltonian 
comprising the remaining terms. The external Hamiltonian H e (see Eq. (6) in the main 
text) in this case is 

H e = -M-{b 2 e 2i " mt - StV 2 *""**) + 2hpftb sin 2u m t, (G3) 
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while the Heisenberg equations of motion read 

_a 2 = - -a 2 - iVL m [b + b\e 2iulmt J + \/Ka 2 , m , 
jb = -zO m (a 2 + a\e 2lWmt ^ + iF(t)e iuJmt + - 2i/%sin 2u m t. (G4) 

We proceed to solve these equations in the Fourier domain using the same techniques 
described in Sec. [F] Specifically, we truncate terms containing frequency components 
to + 2nuo m (integer n) at |n|>2 and solve for a(cu), b(cu) in terms of F(u) and a m (c<j), from 
which any correlation functions for the cavity field or mechanical motion can be obtained. 
The input field is assumed to be in the vacuum state. Similarly, the properties of the output 
field can be obtained from these solutions by using the relation a out = — a m . 

We are specifically interested in the properties of the operator X + out (c<j = 0) = a out (u = 
0) + al ut (co = 0). The general solutions of Eq. ( IG4I) in the Fourier domain are quite cum- 
bersome, so we consider the simplified limit where we set T = n, and take the parametric 
driving strength to be (3 = ^-(1 — St), where 5t<Sl is a small parameter that characterizes 
how far one operates from threshold (J3— >T/2). Expanding to lowest order in K/u m and 5 t 
and ignoring fast-rotating terms, we find the following variance, 

A4, out (. = 0)*44 + + ^(l + S t ) + f . (G5) 

16^4 32^4 k 4 

In particular, at threshold (8t = 0), one recovers Eq. (9) of the main text. Maximum 
squeezing of the variance on threshold is achieved when k = 2(2</>/5) 1//3 u; m , in which case 

(AXi )OUt ) min = (3/2)(50 2 /2) 1/3 . 

Now we consider the effects of cavity loss on the maximum achievable squeezing of output 
light. Starting from Eq. (jH]) for the squeezing at threshold in an ideal cavity (with r~/t), 
we model cavity losses via a beam splitter transformation with the ideal squeezed light and 
vacuum as the two inputs. The output light exhibits reduced squeezing due to mixing with 
the vacuum, given by 

-2 I n\\ ( 1 K '\ ( Ztfr^m , 5 k 2 \ i k' 



(A*U(« = o)W=(i--)(ir + I5;4 + «> < G6 » 



where k', k denote the scattering/absorption loss in the cavity and the total cavity linewidth, 
respectively. In the relevant regime where k'/k<^.1, we can approximate 1 — k'/k~1 
and the squeezing is optimized for the choice k = 2(2/5) 1 / 3 (0 + K'/(2t<; m )) 1 / 3 , for which 
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(AJfi out )m in Ri2.O4(0 + fi;'/(2cc; m )) 2/ ' 3 . We now must choose a set of realistic cavity param- 
eters where this optimized squeezing can be realized, and where r~K is consistent with ( 
being small. As an example, taking a cavity length and waist of L~2 cm and u>~10 /zm, 
k! corresponding to 1 ppm losses per round trip, and sphere parameters r~35 nm and 
dj m /(27r)~0.65 MHz, we find that T~k corresponds to a value £~l/4, which yields squeez- 
ing of ~15 dB in the output light. 

Thus far, we have neglected to consider corrections due to a possibly large position 
uncertainty Ax for the CM motion of the sphere. Specifically, as one approaches threshold, 
one quadrature of motion becomes infinitely unsqueezed, producing a large Ax. At the same 
time, faithful quantum state transfer requires a linear opto-mechanical coupling, where 0(x 2 ) 
shifts in the cavity cooling mode frequency can be ignored. Specifically, the Lamb-Dicke 
parameter 7]=kAx^l for the trapped sphere must remain small. To quantify this effect, we 
consider the situation where we operate away from threshold by an amount that decreases the 
squeezing by just 1 dB relative to (AX^ out ) min . The value of 5 t corresponding to this 1 dB 
increase can be obtained by solving Eq. (1G5j) . and plugged into the solutions of Eq. (1G4I) to 
numerically find Ax. For concreteness, here we associate Ax with the position uncertainty 
in the unsqueezed quadrature of motion. The corresponding Lamb-Dicke parameter as a 
function of sphere size is then plotted in Fig. [6] for the choice cu m /(27r) = 1 MHz, and it is 
seen that i] < 1CF 2 over the entire parameter regime. 
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FIG. 1: a) Illustration of dielectric sphere trapped in optical cavity. The large trapping beam 
intensity provides an optical potential U op t(x) that traps the sphere near an anti-node. A second 
more weakly driven cavity mode with a non-vanishing intensity gradient at the trap center is 
used to cool the motion of the sphere, b) Energy level diagram of mechanical motion (denoted 
m) and cavity cooling mode (ph). The mechanical mode has frequency u m , while the optical 
mode has frequency L02 and linewidth k. Photon recoil induces transitions between mechanical 
states \n m )— >|(n±l) m ) at a rate R n ^ n ±\ (i?o->i shown by dashed gray arrow). The cooling beam, 
with effective opto- mechanical driving amplitude £l m , induces anti-Stokes scattering that cools the 
mechanical motion and allows for quantum state transfer between motion and light. This beam is 
also responsible for weaker, off-resonant heating via Stokes scattering, c) Mechanical frequency u m 
as a function of trapping beam intensity. For all numerical results, we take A = 1 //m, p = 2 g/cm 3 , 
and a high index material (f^~l). d) Optical trap depth Uq (in K) as functions of trapping beam 
intensity and sphere radius r. 
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FIG. 2: a) Mean phonon number (nt) (black curve) versus cavity finesse F (F = itc/2kL) under 
optimized cooling conditions. The sphere size is r = 50 nm and the cavity has a length and waist 
of L = 1 cm, w = 25 /im, respectively. The red curve denotes hf >Ta i n , the fundamental limit of 
cooling imposed by sideband resolution, b) Solid blue curve: optimized EPR variance between two 
levitated spheres, as a function of squeezing parameter e~ 2R . System parameters are identical to 
a). Dashed curve: EPR variance in limit of perfect state transfer, Aepr = e~ 2R . Green curve: 
cavity finesse corresponding to optimal EPR variance, c) Optimized variance (AX^_ out ) m i n (in dB) 
of squeezed output light from an ideal cavity, as a function of sphere size. 



30 



O "~0 0.2 0.4 0.6 0.8 1 




O 0.2 0.4 0.6 0.8 1 




FIG. 3: Comparison of optical forces (in arbitrary units) acting on a dielectric sphere of permittivity 
e = 2 as a function of position kzo- The four figures shown correspond to sphere sizes of p=ky/er = 
0.5,1,1.5,2. The black curve indicates the results obtained from an electrostatic, point-dipole 
approximation, while the red curve denotes exact electrodynamical theory. 
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FIG. 4: Internal temperature of sphere (in K), as functions of background gas pressure and intra- 
cavity intensity. Material parameters for the sphere are given in the text. Optical losses for the 
sphere are assumed to be a) 10 dB/km, b) 100 dB/km, and c) 1000 dB/km. 
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FIG. 5: The number of coherent oscillations N osc allowed before a quantum jump due to shot noise, 
as a function of cavity finesse. The system parameters are given in the text. 
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FIG. 6: The Lamb-Dicke parameter r) = kAx corresponding to the squeezed motional state of the 
sphere, as a function of sphere size. The squeezing parameters are chosen such that the squeezing 
in the output light is increased by 1 dB over (AX^_ out ) m in- The physical parameters of the system 
are taken to be A = 1 /im, p = 2 g/cm 3 , f^~l, and uj m /{2ir) = 1 MHz. 
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